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Abstract 

In a previous paper [Phys. Rev. A 65, 032721 (2002)] we re-examined a model describing the 
structure of the low-energy Ps-H resonances as being due to quasi-bound states of the positron in 
the perturbed Coulomb potential of the H~ ion appearing in the closed, rearranged channel. In 
particular, we wished to understand why the lowest p-state resonance was so far away from the 
lowest quasi-bound (2p) state. We found that the lowest resonance actually corresponds to the 
first excited [3p] state, while the lowest state is not recognizable as a resonance. In the present 
work we repeat our analysis, but this time for the lowest d state. We find that the lowest [3d] state 
does correspond to a resonance shifted moderately. 
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There has been much recent theoretical interest in the positronium-hydrogen [Ps-H] scat- 
tering system [1], and special interest has focused on the properties of elastic scattering 
resonances that, ha,ve been predicted by a variety of theoretical techniques [2, 3]. In par- 
ticular, the model of resonances due to hydrogenic bound states in the re-arranged [e + -H~] 
channel predicts the existence of an infinite sequence of resonances in all angular-momentum 
states [4]. We may expect these resonances to lie fairly close to the levels of hydrogen al- 
though perturbed by the non-Coulomb parts of the negative hydrogen ion’s potential. In 
addition, the L-degeneracy of the hydrogen levels should be broken by this short-range repul- 
sive potential so as to raise the low-L states more than the higher ones, since they are more 
penetrating. A simple S-wave calculation based on this idea was carried out long ago [4], 
including the coupling between the hydrogenic bound states and the scattering continuum. 
The lowest S-state was not expected to correspond to a resonance, since there is a single 
S-wave bound state in Ps-H, but the 2s state did generate a resonance that agreed well 
with other calculations [2]. Based on this success, it was assumed that the correspondence 
between hydrogenic bound states and resonances would carry over trivially to higher values 
of the angular momentum. 

Therefore it was a surprise when the lowest P and D resonances were found [3] to be 
radically shifted from their expected positions, based on the hydrogenic model. That was 
our motivation for re-examining that model, including P-states in our analysis. We were 
again surprised to find [5] that the 2 p state was broadened out of all recognition, and that 
it was actually the 3 p state that produced the lowest P resonance, at a position in good 
agreement with the results of Ref. [3]. Since we now know that the lowest state in both S- and 
P-waves is ineffective at producing scattering resonances, we wish to find out whether this 
situation is true for higher-L waves too. In this Brief Report we outline the corresponding 
calculation for L=2. In this case, we find that the lowest state does in fact correspond to 
the lowest resonance, in contrast to the lower-L cases. 

The method described in Ref. [4, 5] is essentially a two-state approximation: 

y(x,fi,r 2 ) = + {1 -H- 2} + G(x)$(r u r 2 ) (1) 


where 


G(x) = ^-P 2 (x ■ z) 


X 


(2) 
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and the ground states of Ps and H atoms have the following forms: 

e~ r 

The function F(R) represents the Ps-H scattering, both atoms being unexcited, while G(x) 
is the wave function of the positron bound in a d-state to the negative hydrogen ion. The 
system is in the singlet electronic spin state (corresponding to the symmetry of the ionic 
bound state <h(ri, r 2 )) and is explicitly antisymmetrized. 

In Ref. [5] the details of the variational derivation of an integro-differential equation for 
F(R) are described. This has the following form: 

-i[V 2 + k 2 }F(R) + f dR'K(R, R')F(R') + 2 £ = 0 . ( 4 ) 
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Two major simplifications are then made: First, the integral kernel coming from the ex- 
change of the two electrons is replaced by an effective separable form 

K(R, R!) = be~^ R+R ^RR'P 2 {R ■ z)P 2 {R! • z) (5) 


e -pp 

M = ~j==, i>(r) 


and second, the hydrogen ion wave function is replaced by its simplest approximation 

*(ri ,f 2 ) = —e~ z ^\Z=^. 
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The first step is to omit the closed-channel function G(x) which removes the optical 
potential terms from Eq. (4) and solve the resulting static-exchange equation: 


-u"(R) - k 2 u(R) + + Xe~ aR R 2 J dR! R ,2 e~ aR ' u(R') = 0 , 


(7) 


where A = 87r6/5 and F(R) = u(R)P 2 (R • z)/R. This equation has the following solution: 


u(R) = j ( kR ) — tan 77 


n(kR) + e~ aR J2 C n R n 
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where 
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( 8 ) 
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and 


. 3 sin 2 3 cos 2 

j(z) = — 5 sin 2 

z 2 2 

. . 3 cos z 3 sin 2 

n{z) = f- cos 2. 
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The effective range formula takes the form 

k 5 cot rj = — — ■ + \ 4. ^ [3 a 8 + 28a 6 k 2 + 178a 4 k 4 + 28a 2 k 6 + Sk 8 ] . (10) 

By fitting this expression to the L — 2 phase shifts calculated in the static-exchange approx- 
imation Ref. [6] we find a best fit for a = 1.30244 and A = —1.02029, and we will assume 
that the separable form of the exchange kernel shown in Eq.(5) is appropriate even when 
we add the correlated terms coming from the closed-channel hydrogen-ion parts of the wave 
function of Eq.(l). 

The virtual bound state wave function G(x) can consist of one or more D-state eigen- 
functions of the differential equation 

6 

-9Ndi x ) + — 9Nd{x) + U(x)g Nd {x) = e Nd g Nd (x), (11) 

X 

where U(x) is the potential due to the simple form of $ shown in Eq.(6): 

v{x) = -l + ie -™*(i + Z ). (12) 


The two lowest eigenvalues are e 3d — —0.1097 Ry and — —0.061 8Ry, which lie very close 

to the unperturbed hydrogenic values of —1/9 and —1/16. We are interested in finding how 
close to these positions the corresponding scattering resonances lie. (Note that the second 
of these lies above the inelastic threshold of the Ps atom at e = —3/8 — Eh- + Eh + Ep s = 
-0.0695.) 

We follow the same procedure in calculating the separable potentials V Nd {R) as for the 
earlier S- and P-wave calculations of Ref. [5]: 


V Nd (R 1 ) = J J dp l df^{p l )ii{r 2 ) ^ N -^ P 2 (x-z)^r l ,r2) 


. 2 , . 2 2 

&Nd H U(x) 

x Pi P2 


(13) 

Here A = &h- + e Nd — Eh — Ep s — k 2 / 2, so that for the lowest resonance A 3d = 0.3348 — 
k 2 / 2. Once the potentials have been evaluated the numerical solution of the scattering 
equation by the noniterative method is quite straightforward. 

The potential Vpr d (Ri) is the result of an integration over the two remaining independent 
Jacobi coordinates describing the system, pi and f 2 . But there are additional coordinates 
involved in the integration, especially the argument of the Legendre polynomial giving the 
angular dependence of the bound state function G(x). We expect the angular dependence 
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of the optical potential to involve only D-states of the variable pr = Ri • z, so the variables 
in the integrand must be appropriately transformed. (The integration over df 2 is trivial.) 
Since x = R\ -j- f>i/2 we can write 


X ‘ Z = ^ (R\PR+^PiP p 


X 


(14) 


where we define p p = pi ■ z and p ~ pi ■ R\- The first term in Eq.(14) has the desired form, 
and we can write the second in terms of 

A/jr _j_ j 

Pp = -r £ Y^pr^Y-^ 4>) = MkA - V (1 - /4)(1 - M 2 ) cos(^ - <j)) (15) 

13 m =- 1 

In our previous work on P-states, at this point we could neglect the ^-dependent terms 
which would vanish after integration over d(f>. But in the present case we must construct 
P 2 (x ■ z) which is quadratic in its argument, and the quadratic ^dependent terms do not 
vanish after integration. After some algebra we obtain the transformed expression 




x z 


pi 


R x + Ripift + —Pzift) 


(16) 


Finally, we can now write the potential explicitly in terms of a double integral: 

-Zt9{x) 


V Nd {R) = 4^ {— P 2 (p R ) j dpp 2 e~* J dpe~ ZT ^Q (A Nd + P). (17) 

o -1 


Here we have used the definitions 
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Q = R 2 + Rpv + -p 2 P 2 {p) 


which comes obviously from Eq.(16), and 

P = e~ (1 + Z)x (J + 1 + Z^- U(x) - 

Everything has now been re-written in terms of the three variables R, p, and p, keeping 
in mind that [x, r] = \J R 2 + fP/A ± Rpp. In practice we separate V Nd (R) into two parts, 
one independent of energy and one proportional to A Wd , carry out the double integrals 
numerically for a range of values of R , and then fit the results smoothly by interpolation. 
The resulting potentials are then used in the scattering equation (4), after expansion in 
partial waves, to compute the D-wave phase shift as a function of energy. 

The results are not as surprising as were those presented in Ref. [5]. In Fig.l we show the 
phase shift when only the 3 d closed-channel state is included in G(x). In this case a single, 




FIG. 1: Ps-H D-wave phase shifts. The squares are the result when only the 3d state is included, 
and the crosses are the result when the 4d state is added. The curves are best-fit analytical 
functions of the standard form [Eq.(18)]. The two horizontal lines are at tx and 2i r, respectively, to 
guide the eye. 


fairly narrow resonant feature is seen, and when it is fitted with the standard resonance 
formula 


r](E) — A + BE + arctan 


(18) 


.2 (E r - E) 

we find the results shown in Table I. (The energies there are the Ps-H scattering energies 
k 2 /2, not energies measured relative to H _ ). You can see that the resonance energy obtained 
this way is very close to the result in [4], although the width is too large. 

Also in Fig.l we plot the phase shift when we add the Ad state. Here we find two 
resonances, which lie quite close together. The results of a two-term fit are also shown in 
Table I, although the higher one is perhaps not too significant since it lies above the Ps 


6 




TABLE I: D-wave resonances in various approximations. The energies e^d do not include coupling 


to tne continuum. 


x lie lowest lesonance is labeled 1 and the second is 2. 


Method 

Ps scattering energy (eV) 

Width (eV) 

e 3d 

4.555 


e 4d 

5.207 


Ei(l term) 

4.708 

0.394 

Ej (2 term) 

4.729 

0.327 

Ei [4] 

4. 710T0.0027 

0.09 25T0.0054 

E 2(2 term) 

5.253 

0.266 


inelastic threshold of 0.375 Ry and has not been previously calculated. With the inclusion 
of this second state the position of the first resonance rises slightly and is in somewhat poorer 
agreement with Ref. [4], while its narrower width is a slight improvement. 

We conclude from these calculations that the D-wave resonances behave almost exactly 
as we would have expected before the problem with the P-wave scattering arose. That is, 
coupling with the continuum shifts the position of the lowest resonance, but only by a modest 
amount. There is no sign of the phenomenon that startled us in the P-wave work, where 
the lowest potential resonance disappeared due to that coupling. We now feel confident 
that calculations of higher-L resonances (which we do not intend to carry out) will continue 
to act normally, although their positions above the inelastic threshold would necessitate 
multi-channel formulations. 
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